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Abstract 

We compute the thermal photon spectrum in the Pb+Pb collisions at the CERN 
SPS energy using thermal emission rates and a hydrodynamic description for the evo- 
lution of produced hot matter and compare our results with the measurements of the 
excess photons by the WA98 collaboration. Our results show that the measured photon 
spectrum can be reproduced with realistic initial conditions which take properly into 
account also the finite longitudinal size of the initial collision zone and which simul- 
taneously describe well both the transverse and longitudinal hadron spectra. In the 
scenario with initial formation of QGP the recently calculated emission rate, complete 
to order a s , reproduces the measured spectrum. However, the experimental spectrum 
can also be reproduced in a purely hadronic scenario without transition to QGP state, 
but a high initial temperature, much over the values predicted for the phase transition 
temperature T c , is required. 
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1 Introduction 



For the study of the properties of matter created in a high energy heavy-ion collision, 
photons and lepton pairs provide the most direct probe since, after being emitted, 
they escape the system without further interactions. Unfortunately, electromagnetic 
interactions occur during all the different stages of the collision, ranging from those 
between the incoming partons to the decays of final hadrons, and it is difficult to 
disentangle experimentally the contributions from these different sources. 

After a painstaking analysis, the WA98 collaboration has succeeded in isolating an 
excess of photons over those originating from the decays of final hadrons [l], [| . Since 
a QCD calculation of photon production from the primary interactions of incoming 
partons fails to reproduce this excess^] [0, |3|, they can be associated with the secondary 
interactions of final state particles. In a phenomenological analysis of heavy-ion colli- 
sions, based on the assumption of the formation of a thermal system, these are often 
called the thermal photons. 

Since the publication of the WA98 data there has been several attempts to describe 
the data using a thermal description for the expanding matter || || |7|, |[ |9j and the 
calculated thermal emission rates of photons. In these studies the time evolution of 
matter has either been given by a simple parameterization of spherical expansion Q 
or described by using boost invariant hydrodynamics |6|, [7], |J. In the papers ap- 
plying hydrodynamics, the reproduction of the data has required either large initial 
temperature which led to unrealistically short initial time for boost invariant flow ||, 
significant initial radial velocity |6], [7], |8|], or large modifications in hadron masses |J. 

In the papers applying hydrodynamics || || [Tj, || the validity of boost invariant 
expansion is taken at face value. However, from the hadron data at the SPS, we know 
that the rapidity window where the expansion is boost invariant is narrow if it exists at 
all. Even if the flow is approximately boost invariant at later times, it can be expected 
to deviate from the scaling behaviour for r<l fm/c because of the finite longitudinal 
extent of the collision zone as discussed in more detail below. It is also known ]ITJ 
that if the assumption of boost invariance is relaxed and the longitudinal expansion is 
constrained to reproduce the measured rapidity spectra of hadrons, the initial energy 
densities become strongly peaked and larger by a factor two than the Bjorken estimate. 
As is evident in || the thermal photon production at large is very sensitive to initial 
temperature and dominated by the earliest moments of the expansion. Therefore, the 
photon production even at midrapidity becomes sensitive to the details of the initial 
state including its longitudinal structure. 

In this letter we take a conservative approach to photon production. We use well 
known vacuum properties of mesons also in dense medium and set the initial radial ve- 
locity to zero. Instead of studying these effects we concentrate on the effects on photon 

1 It should be kept in mind that a consistent description of inclusive photon data in hadronic 
collisions in fixed target experiments is difficult to achieve 0. 



1 



spectra caused by different initial states and equations of state (EoS) of the expanding 
matter. In all cases the calculation is constrained to provide a reasonable description 
of both transverse and longitudinal hadronic spectra JTTJ, not only the final particle 
multiplicity as in 0, || |8j . In our model finite baryon number density is included in the 
EoS. We assume local chemical equilibrium with zero strangeness density until kinetic 
freeze-out. This allows a reasonable reproduction of charged particle and net proton 
spectra but fails to describe the details of strange particle yields. However, this is not 
a severe problem since reproduction of negative particle spectra means the description 
of the bulk of the matter is realistic. The processes involving strange particles are not 
included in the photon production rates in hadron gas either. To compare with the re- 
sults in literature, we calculate the photon spectra also with boost invariant expansion. 
As thermal production rates we use rates provided in the literature. For the QGP we 
use the recently obtained full order of a s result [0, [13| for zero quark chemical poten- 
tial, n q = 0. This result contains both the lowest order one-loop rates 0, [15| with 
exact non-logarithmic term and completes the earlier two-loop calculations 
through summation of leading contributions from multi-loop terms up to all orders 
hadron gas we use the parameterizations 
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In 

IB), HH °f the rates calculated using effective 



Lagrangians and vector-dominance model (14], |19|. 



2 The framework for calculating thermal spectra 



Consider first the emission rates of photons from a thermal matter. Since we assume 
that the photons escape the system after emission the emitted spectrum is not the 
thermal, black-body radiation spectrum. 

In the QGP one is faced with the problem of infrared singularities in the rate cal- 
culation. In lowest order this leads to a logarithmic energy dependence on uj/T in 
addition to the usual exponential (in Boltzmann approximation) exp{— cj/T} depen- 
dence |L4], [15| . The next-to- leading order 2-loop result [TJ|, [T| turned out to be of 
the same order in a$ and of the same size or larger than the leading 1-loop result. It 
was argued that also the leading parts of the higher-loop contributions are of order 
a s [20]. Recently the resummation of order a s contributions from all orders in the 
multi-loop expansion was achieved, fully completing the order of a s analysis of the 



photon emission rate in the quark-gluon plasma [JT2|, [13 
the parameterization of the rate provided in [JT3 



In our calculation we use 
It should be remembered that this 



result is for zero baryochemical potential whereas in the collisions at SPS energy net 
baryon density is relatively high. In the leading order, the effect of non-zero fi q has 
been considered plj and, in good approximation, it was found to be included into the 
rate by replacing the multiplicative factor T 2 by T 2 + fig/ir 2 . We have checked that 
in the plasma phase the ratio of these factors is typically 1.02 < 1 + /i 2 /(7rT) 2 < 1.05 
indicating only a few per cent increase to the calculated plasma contributions. 
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Another problem in applying the calculated rates in nuclear collisions is that the 
assumption of high temperature limit with an ideal QGP may not be valid. The lattice 
QCD calculations of energy density and pressure indicate that the ideal gas or Stefan- 



Boltzmann limit is approached slowly with increasing temperature [22|. On the other 



hand recent lattice calculations on lepton pair emission rate |23fl show that, at least in 
some parts of the phase space, the lowest order (in a s ) perturbative result is quite close 
to the full lattice result. This is encouraging for the type of phenomeno logical study 
as presented here and supports the interpretation of data in terms of results based on 
perturbative rates. 

In the hadron gas the scattering processes tctt — > p7 and up — > 7r7 are included in 
the photon emission rate. In addition to the contributions from interactions described 
with a pseudoscalar- vector Lagrangian fifM, the process np — > ivy gets a large contri- 



bution from the irp interaction with the a± axial meson: irp — > a\ — > 7T7 In our 

calculations we use the parameterizations provided by the authors in |18|, [19[ for these 
rates. Since in the WA98 analysis the decay contributions of final state hadrons are 
subtracted, only those p — > 7r + 7r + 7 and uj —>■ 7r° + 7 decays which occur in the 
hadron gas before freeze-out are included. Even though the decays of final 7r°'s and 
77's dominate the decay photons from final hadrons, their decays are so much slower 
than those of p and u that they can be neglected during the hadron gas phase. The 
effects of baryons on the emission rates are not known and our results contain rates 
from mesonic processes only. In the hydrodynamic calculation of expansion baryons 
are included. 

We study the evolution of the system using three different equations of state (EoS). 
Two of the EoSs describe a first order phase transition from hadron gas to QGP with 
transition temperatures 165 MeV (EoS A) and 200 MeV (EoS D) and the third one 
describes the matter as a hadron gas also at high temperatures (EoSH). The QGP 
is an ideal gas of quarks with 2.5 flavours and gluons. The hadron gas contains all 
known hadrons with masses up to 2 GeV. A more detailed discussion can be found 
in [[11], |24| . The present results from lattice calculations for the transition temperature 
are T c (n f = 3) = 154 ± 8, T c (n f = 2) = 173 ± 8, and T c (n f = 2) = 171 ± 4 MeV, 
where the first two values are for improved staggered fermions and the last one for 



Wilson fermions ||22|| . For this reason we consider the EoS A with T c = 165 MeV as our 
reference EoS. 

As described in the introduction we want to study the effects of initial conditions 
on the longitudinal expansion and the photon production. For that purpose we do the 
calculation using three different descriptions of initial state. For two of them, labeled 
IS 1 and IS 2, the longitudinal expansion is not boost invariant. To be able to distin- 
guish what is the effect of longitudinal expansion, we also carry out the calculations 
using the boost invariant expansion scenario (labeled BI in the following). When the 
expansion is not boost invariant, we assume that the rapidity distribution of energy in 
the initial state is Gaussian. The width of the rapidity distribution increases with the 
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IS1 


IS 2 


BI (0.3) 


BI(l.O) 




EoSA 


EoSH 


EoSA 


EoSH 


EoSA 


EoSA 


T max (MeV) 


325 


275 


265 


245 


364 


244 


T(z = 0) (MeV) 


255 


234 


214 


213 


301 


214 



Table 1: Maximum value and average initial temperature for different EoS and initial 
state. Average temperature is calculated at z — 0. 



transverse coordinate r since the nuclei become more transparent towards the edges. 
The normalization is fixed by requiring the energy per unit transverse area in the ini- 
tial state, e(to,r) = JdzT 00 (to,r, z), to equal the energy per unit transverse area of 
incoming nucleons JTT| . 

To convert this distribution to spatial energy density distribution we have to spec- 
ify the z dependence of the initial longitudinal velocity. This choice is by no means 
unambiguous and since the choice of the velocity profile affects the maximum initial 
temperature and thus the photon yield, we carry out the calculations using two dif- 
ferent profiles. In both cases the longitudinal size of the system is constrained to 3.2 
fm and the longitudinal flow velocity at the edge of the system is similar. In the case 
of IS 1 we assume that the longitudinal flow rapidity increases linearly with distance z 
from the center of the system. In our second choice IS 2 the dependence is nonlinear 
and the velocity is smaller at small z, but becomes at large z larger than in the linear 
case IS 1 (see ref. |Tl|). In the case of IS 1 the energy density is strongly peaked in 
the middle whereas IS 2 leads to an almost flat distribution. Consequently IS 1 leads 
to a larger value of the maximum initial temperature than IS 2 (see table [I]). I n the 
following we use IS 1 as our reference expansion scenario. 

In the case of boost invariant longitudinal flow, v z — z/t, we assume the transverse 
distribution of energy density e(r ,r) to be proportional to the number of participants 
per unit transverse area. Its normalization is chosen to reproduce the observed pion pt 



distribution in most central collisions |Pq| . The normalization of e obviously depends 
on the initial time To which is the time for all the matter in the given rapidity range 
to reach local thermal equilibrium. In the boost invariant scenario this represents the 
main uncertainty of the calculation. It should be observed that the hadron spectra are 
not very sensitive to To as long as tq is much shorter than the transverse time scale 
VSRa and if the initial densities are chosen to give the same total entropy. However, 
since the photon rates have a strong temperature dependence, especially the high-/c T 
part of the thermal photon spectrum changes considerably with different choices of To. 
To illustrate this effect, we carry out the boost invariant calculation for two different 
initial times, r = 1.0 and 0.3 fm/c, labeled BI (1.0) and BI (0.3), respectively. In both 
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cases we fix the initial densities by taking the density profile used in |25j at tq = 0.8 
fm/c and scale it to initial times r = 0.3 and 1.0 fm/c by requiring that entropy per 
unit rapidity does not change. 

At finite collision energy the longitudinal size of the collision zone is ~ 2R J ^/ / y cm which 
is also the measure of the time interval for the primary particle production. Since in 
the boost invariant expansion the longitudinal size of the system is ~ 2r, the initial 
time should satisfy To>-Ra/ 'c7cm ~ 1 fm/c for a consistent use of scaling hydrodynamics 
at SPS.Q Even though we do not consider initial times as short as r = 0.3 fm to be 
realistic for a boost invariant expansion at the SPS, they may be useful in estimating 
the limits of photon spectra. In scenarios without boost invariance, the essential initial 
parameter is the longitudinal size of the system. In fixing the longitudinal size from the 
longitudinal collision geometry we have implicitly assumed that local thermalization 
times are shorter than this longitudinal scale. 



3 Results 

We have calculated the photon emission for all the combinations of three equations 
of state, EoSA, D and H, and the initial conditions, IS 1 and IS 2. For comparison 
we have performed the calculation also for the initial conditions BI(0.3) and BI(l.O) 
with EoS A. To illustrate the quality of our hydrodynamic description of the nuclear 
collision, we show in fig. [I] the measured ir° spectrum with our calculated results for 
some of the cases we have studied. The solid line (IS 1, EoS H) and the dashed line (IS 2, 
EoS A) are the two extreme cases for the calculations with the non-scaling longitudinal 
expansion. All other combinations of non-scaling flow and EoS give spectra within these 
limits. The dotted (BI (0.3), EoS A) and dashed-dotted (BI (1.0), EoS A) lines represent 
calculations with boost invariant longitudinal flow, EoSA and initial conditions from 
25j scaled to To = 0.3 and 1.0 fm, respectively. We have not tried to fit the 7r° spectrum 
but, instead, used the same initial conditions as in earlier studies of hadron spectra 
from several SPS lead- lead experiments |JTTJ, 

In fig. 2 we show the contribution from the hadron gas (dotted line) and contributions 
from plasma using the 1-loop rate (dashed-dotted line), the sum of 1-loop and the order 
of a s part of 2-loop rates (dashed line) and the total order of a s rate obtained summing 
the contributions from multi-loop calculations up to all orders. The initial conditions 
IS 1 is used in the calculation. We conclude that the emission from plasma is dominant 
for the large transverse momentum range, > 2 GeV. The suppression of the 2-loop 
result due to the resummation of higher loop contributions is roughly 30 % and the 



2 In p6[ a saturation criterion for minijet production is presented leading to a time scale ~ 0.3 fm/c 
at the SPS. This is the time scale for the production of a minijet of transverse momentum ~ 0.7 GeV 
and, at the SPS, cannot be taken as the time interval for the formation of all minijets. Indeed, it is 
not clear if the saturation applies or if soft processes with a longer time scale are more significant at 
the SPS. 
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Figure 1: Calculated it ^-distribution for 
different initial states and EoSs compared 
with the measurement by the WA98 col- 
laboration For details, see the text. 



Figure 2: Thermal photon emission (IS I 
and EoS A) compared with the WA98 
data ffl. The contribution from plasma 
is shown for rates calculated in leading 
order (1-loop), next to leading order (2- 
loop) and in order a s with contributions 
from multi-loop calculations in all orders 
(total). Photons from hadronic phase are 
shown separately (HG). 



1-loop contribution is typically ~ 40 % of the total QGP rate. The shapes of the three 
plasma curves are remarkably similar. 

In fig. 3 the total photon spectra are depicted for different initial conditions using 
EoS A in each case. As discussed in chapter 2, the initial values of energy density or 
equivalently of temperature are smaller for initial conditions IS 2 than for IS 1. In case 
of IS 2 the photon emission is reduced as compared to IS 1, especially at larger values 
of &>r- At later times the evolution is quite similar in both cases and the contributions 
from hadron gas do not differ significantly. For IS 2 the high-fc^ end of the spectrum 
is not reproduced. 

We show in fig. 3 also the results from calculations with longitudinal boost invariance 
starting the evolution either at r = 1 or 0.3 fm/c with the same entropy per unit 
rapidity, dS/dy, in both cases. With this choice of initial times the results follow 
closely the nonscaling calculations so that for tq = 0.3 fm/c (BI (0.3)) the spectrum is 
close to that from initial conditions IS 1 and for tq = 1 fm/c (BI (1.0)) to that from IS 2. 
The high-fc^ data is reproduced only with high temperature initial states, BI (0.3) for 
the longitudinal scaling flow and IS 1 for the non-scaling case. Since the cooling due to 
the longitudinal scaling expansion is very fast in the case of short initial time, a much 
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Figure 3: Thermal photon emission for 
different initial states. Data are from |l[ 
and calculations are done using EoS A and 
total order of a s rates in the plasma. 




k T (GeV/c) 



Figure 4: Thermal photon emission for 
different EoS and initial state with the 
contribution of initial pQCD photons (up- 
per set) and without (lower set). The data 
and results of the upper set have been 
scaled by a factor 100. 



higher initial temperature is needed in case of BI (0.3) than in case of IS 1, see table 0. 

As was argued in the previous chapter, we do not consider the use of a short initial 
time, such as r = 0.3 fm/c, to provide a fully consistent description of the early 
stage of a heavy-ion collision at the SPS energy even though the flow at mid rapidities 
at later times might evolve close to scaling flow. We have seen, however, that the 
early hot stage is essential for photon emission emphasizing the need for its reliable 
description. One could argue that the longitudinal size of the collision region is reduced 
below ~ 2 fm by shock compression but this would also slow down the initial flow from 
the scaling velocity v z = z/t. In fact, this would lead to an initial state with even 
higher temperatures and slower expansion velocity than the one described by IS 1 and 
consequently to results overshooting the measured photon spectrum. 

Part of the measured prompt photon spectrum is expected to come from the primary 
parton-parton interactions of incoming nuclei. A perturbative QCD analysis, which 
includes also a study of the effect of intrinsic transverse momentum of partons, gives 
results which are not more than 30 % of the measured spectrum when (p^)mtr — 0.9 
GeV 2 ||. If no intrinsic partonic pt is included, calculations are an order of magni- 
tude below the data. A comprehensive analysis of prompt photons in hadron-hadron 
collisions at comparable energies does not favour large values of intrinsic (pj) |4j]. It 
looks unlikely that the WA98 photon data can be explained in terms of pQCD photons 
from interactions of incoming partons but it would be premature to claim that this 
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possibility is definitely excluded. 

In fig. |4] we show the photon spectra for EoS A using initial conditions IS 1 and for 
EoSH using both initial conditions IS 1 and IS 2. We have performed the calculations 
also for the EoS D with T c = 200 MeV. Even though the different contributions change, 
the total spectra are qualitatively similar for EoS A and EoS D|. The reason to this is 
that at high the highest temperatures dominate and they are not changed. The 
higher production rate of photons in hadron gas than in plasma Jl7| leads to a higher 
yield for EoS H and IS 1 than for EoS A and IS 1 even though the initial temperature 
is higher in the latter case (cf. table p. The spectrum reaches the upper limits of 
experimental data points around kr ~ 2 GeV when EoS H and IS 1 are used. The 
choice of initially cooler initial state IS 2 with EoS H decreases the yield sufficiently in 
kx ~ 2 GeV region but also creates a distinctive convex shape of the spectrum. 

To get a better understanding of the total photon yield in different cases we have 
parameterized the pQCD photon spectrum of || as presented in and added it to 
the thermal yield. Resulting spectra are displayed as the upper set of curves in fig. 4 
for kr > 1.2 GeV. The differences are small both in normalization and shape. We feel 
that with the present accuracy of the data a x 2 fitting would not be useful in ruling 
out any of the curves. The hadronic alternative, EoS H, with high temperature initial 
state IS 1 tends to stay above the data. The hadronic EoS with the lower density initial 
state, IS 2, goes through the data but there is a tendency for the curve to be slightly 
steeper than the data. The spectrum obtained by assuming a phase transition, EoS A 
and initial state IS 1 gives a very good description of data and seems, in particular in 
the shape, to be closest to the data. 

We should like to point out that even with the initial state IS 2 the maximum value of 
the initial temperature in the hadron gas alternative is 245 MeV (cf. table [l]). We have 
checked that for the high-fc^ region, contributions from the short time interval with 
highest temperatures dominate. In all scenarios the contribution from matter with 
temperature T < 200 MeV is roughly 40 % of the total thermal yield at k^ ~ 3 GeV. 



4 Discussion and conclusions 

We have studied in detail the dependence of thermal photon emission in central lead- 
lead collisions at CERN SPS describing hydrodynamically the expansion of produced 
matter and using known emission rates in hot matter. We have considered emission 
both from QGP and hadron gas and compared with each other the different contri- 
butions to the thermal rate. In the hydrodynamic calculations we have used hadronic 
spectra to constrain the initial conditions including the longitudinal structure. The 
significance of initial conditions is emphasized by the fact that most of the photons, in 

3 In rcf. the difference between the yields for EoS A and D was larger than here because we 
used only lowest order production rates. 
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particular at high kr, are emitted in a short time interval right in the beginning of the 
expansion. For comparison we performed also the boost invariant calculation. 

We find that the measured photons can be described in terms of a thermal model 
with hydrodynamical expansion either by using EoS A with phase transition or EoS H 
with hadron gas phase only. The model is constrained to reproduce the transverse 
and longitudinal hadron spectra. However, in describing the photon emission a high 
temperature initial state must be assumed. With the longitudinal scaling expansion 
this is achieved by choosing a short initial time r = 0.3 fm/c which is not in complete 
accordance with longitudinal geometry. The need for high temperature initial state 
can be avoided by assuming a non-zero initial transverse velocity in case BI(l.O). We 
have tried this but, e.g. the use of a linearly increasing initial velocity as in ref. |7j with 
v r (r = Ra) = 0.3, which produces very closely the same photon spectrum as IS 1 or 
BI(0.3) in fig. 3, leads to a n° spectrum well above the case IS 1, EoSH in fig. 1. We 
cannot lower the 7r° spectrum by simply decreasing the initial temperature since we 
would then go below the hadron multiplicities. 

The shape of the photon spectrum clearly depends on whether phase transition is 
assumed or not. Unfortunately this dependence can not be used yet to differentiate 
between scenarios since the differences become as small or smaller than the experimen- 
tal errors when the initial pQCD photons are added to the spectrum. One may argue 
that after the addition of initial pQCD contribution the spectrum with EoS A and IS 1 
still reproduces the shape best, but the differences are small. 

With QGP formation we find the higher-loop contributions essential for reproducing 
the data over the whole observed range. In particular, the shape of the spectrum 
follows the shape of the data and the normalization leaves room for the initial pQCD 
photons which are estimated to be of order 1/3 of the data with spectral shape very 
similar to that of the data. On the other hand, the experimental data rules out, within 
the thermal model, an increase of the rate by a large factor. This could be interpreted 
as lending support for the assumption that the emission rate of photons from plasma 
can be calculated perturbatively already in the temperature range 200 . . . 250 MeV and 
the main part of the rate is given by the full order of a s result W2L IT3 . 



In summary, the WA98 photon data supports a realistic thermal model with high 
density and temperature initial state constructed to describe the main features of 
hadron data. The spectra obtained by assuming a phase transition, EoS A, or pure 
hadron gas, EoS H, deviate somewhat in properties but the accuracy of the present 
data and the present understanding of initial pQCD photons is not enough to rule out 
either alternative. 
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